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Abstract 

The occurrence of chaos for test particles moving in a Taub-NUT space- 
time with a dipolar halo perturbation is studied using Poincare sections. We 
find that the NUT parameter (magnetic mass) attenuates the presence of 
chaos. 
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The Taub-NUT (Newman, Tamburino and Unti) spacetime [|T],0 is one of the most bizarre 
solutions of the vacuum Einstein equations. Because of its many unusual properties it has 
been described as a "counterexample to almost anything" 0. It has closed timelike curves, 
is nonsingular in a meaningful mathematical sense but is not geodesically complete, etc. For 
certain range of the coordinates, it can be seen as a Schwarzschild monopole endowed with 
a "magnetic mass" 

The Euclidean version of this metric has recently received some attention due to the 
fact that it is closely related to the dynamics of two non-relativistic Bogomol'nyi-Prasad- 
Sommerfield (BPS) monopoles [^. The asymptotic motion of monopoles corresponds to 
geodesic motion in Euclidean Taub-NUT space, this motion is integrable. This fact has 
motivated the study of geodesies in Euclidean Taub-NUT and related spaces 0. 

Examples of chaotic motion in General Relativity are the geodesic motion of a test 
particle moving in the geometry associated with : a) Fixed two body problem 0], b) A 
monopolar center of attraction surrounded by a dipolar halo [^,^ (in Newtonian theory 
this system is integrable), c) A monopolar center of attraction surrounded by a quadrupole 
plus octupole halo []TU|, d) Multi-Curzon and multi-Zipoi-Vorhees solutions |]rT|, and e) A 



rotating black hole (Kerr geometry) with a dipolar halo [^. Also gravitational waves can 



produce irregular motion of test particles orbiting around a static black hole |112| , [13 . 

In this Letter we consider the geodesic motion of test particles moving in a Taub-NUT 
spacetime perturbed by a distant distribution of matter that can be represented by a dipole. 
The case of a center of attraction (without NUT parameter) perturbed by a dipolar halo 
was studied in 0, the combined relativistic effects and the breakdown of the reflection 
symmetry in this case produces a non integrable motion. Our main goal in this letter is to 
study the effect that the magnetic mass has on the chaotic motion of test particles. The 



geodesic motion in a pure Taub-NUT spacetime was studied in JT^ , this case is completely 
integrable. Also generalizations and perturbations of this spacetime has been considered 
from the view point of spacetime dynamics 



The metric that represents the superposition of a Taub-NUT metric and a dipole along 



the 2;-axis is a stationary axially symmetric spacetime. The vacuum Einstein equations 
for this class of spacetimes is an integrable system of equations that is closely related to 
the principal sigma model ||T6[. Techniques to actually find the solutions are Backlund 
transformations and the inverse scattering method, also a third method constructed with 
elements of the previous two is the "vesture method" , all these methods are closely related 
||16| . The general metric that represents the nonlinear superposition of a Kerr-NUT solution 
with a Weyl solution, in particular, with a multipolar expansion can be found by using the 
"inverse scattering method" For the particular case of a Taub-NUT metric with a 



dipolar halo we find 

ds^ = gtt{r, z)dt^ + 2gt^{r, z)dtd(j) + g^4>{r, z)d(j)'^ + f{r, z){dz^ + dr^), (1) 

where 

g,^ = _e-2^™[e-^^"(e-«^^(m + ifiu'' - 1) + e-«^(m - ifiu'' - 1) 
^2e-4(-+i)^(n2 - v^)b^) + e-4(-+i)^(e-8^" + l){v^ - l)b^]/H 

g^^ = -[2e-2^("+"+i)(e"^''"(e~'''''(m + l)(n - v){u + 1) 

-e-^^(m - + v)iu- l)){v + 1) + e-^^"(m + l){u + v) 
{u + l){v - 1) - e-^^{m - l){u - v){u - l)(f - l))b]/H 

^2e-4^(-+i)(n2 - v^){u^ - l)b^){{v^ - 1) 
^^-4D{v+i)^^-sDu^^ + 1)4 + {v- iy){u^ - l)b^)/H 

2 

/ = ^Hexp[D^{uh^ -u'^ -v^ + l) + 2D{uv + 2u + 2v + 2)] 
H = e-4^"[e-«^^(m + l)\u + 1)' + e-«^(m - ifiu - if 

= + (2) 

The coordinates (r, 0, z) are dimensionless and have the range of the usual cylindrical 
coordinates. They are related to u and v by: z = uv and r = (m^ — 1)^/^(1 — v'^Y^'^, u > 1 



and —l<v<—l. Our units are such that c = G = 1; m D and b represent the mass, 
the halo dipole strength, and b the NUT parameter, respectively. The constant a represents 
the geometric "sum" of the "electric" and "magnetic" mass. The coordinate transformation 
t' = t, u = R/m — 1, V = cos^, ip' = ip reduces (0) with D = to the Taub-NUT solution 
in the usual spherical coordinates . 

To study the small dipole perturbation case is better to use the metric obtained by 
keeping the first order terms in the dipolar strength D in the exact metric (^) . This approx- 
imation, for the parameters and range of coordinates used, will not produce a significant 
information loss; we shall comeback to this point later. We find for (7^1, = g^^^^^ + Dgjj^^^, 

g^^ = -{2{4m\ - In? + 2muv + u^v - 3t; + 2)Du - F)(l - 

gt^ = 2[-2{{{3v^ - + {2v -3){v + - + v)m - 2[u^ - v^)m^u 
+ {5v^ -2v- l)u^ - {3v - 2) uv — u^)D — F{u^ — l)v\b/ F"^ 

g^^ = -a^[2{2{2{{2v^ - v + 1) ~ {v + l)u'^)m^u + {u'^v + 3u^ + 6u% - 2u^ 

+17 -\){v- l)m - {hv - V)u + - + 4v?v^)F{v + 1) - (2((2w + \)v? + w - l)m 
+ + 2)u^ - 2m^u + huv){'4{u^ + \){v^ - \)mu - 4{u^ - v^)ir? + 2(3t;2 - \)u^ 
+ - \)u^ - - 1))(-1)Z^ + (4(^2 + - \)mu - 4(^2 _ v^)m^ 
+2{3v^ - ly + {v^ - - 3v^ - 1)F]/F^ 

f = a^{2{2{u^ -v + l)m- {3v - 2)u + 2m^u + u^v)D + F) 

F = 1 + 2mu + (3) 

The geodesic equations for the metric (|^) can be cast as 

i = g''E,, <P = g^'E,, (4) 
r = -^[gfEaE, + fAr' - z') + 2/,,ri], (5) 

z = -jj[g-^E,E, + Uz' - t') + 2f,.rzl (6) 

where the dots denote derivation with respect to s and the indices a and b take the values 
{t,(f)), g""^ stands for the inverse of gab- Et = —E and E^ = L are integration constants; 



E and L are the test particle energy and angular momentum, respectively. The set (^-(^ 
admits a third integration constant 

Es = g'''EaE, + f{r' + z') = -l. (7) 

Thus to have complete integrability we need one more independent constant of integration. 
In the case of pure Taub-NUT solution {D = 0) we have a fourth constant due to the 
existence of a third Killing vector and a fifth (Runge-Lenz vector) related to the existence 
of a Killing- Yano tensor. 

The system d^)-® can be written as a four dimensional dynamical system in the vari- 
ables (r, z, Pj. = "T, Pz = z). A convenient method to study qualitative aspects of this system 
is to compute the Poincare sections through the plane z = 0. The intersection of the orbits 
with this plane will be studied in some detail for bounded motions. We shall numerically 
solve the system (|)-® and use the integral to control the accumulated error along the 
integration; we shall return to this point later. 

The Poincare section for different initial conditions with energy E = 0.967 and angular 
momentum L = 3.75 moving in an exact Taub-NUT geometry {D = 0) with NUT parameter 
b = 0.28 and "total mass" a = 1 are presented in Fig. 1. We have the typical section of an 
integrable motion, i.e., the sectioning of invariant tori, for integrability and KAM theory, 
see for instance [|19|. The values of (7,E and L will be kept unchanged in our numerical 
analysis. 

The motion of test particles around a Kerr and a Schwarzschild black hole with a dipolar 
halo (6 = and D ^ ) is chaotic and it was studied in some detail in for different 
energy shells. In Fig. 2 we show the Poincare section for D = 0.0005 and the same values of 
E = 0.967 and L = 3.75 as in Fig. 1. We find islands of integrability surrounded by chaotic 
motion. The two isolated islands around the points (10, 0.05) and (5, 0.075) correspond to 
the same torus. In the case studied in 0] they were closer. 

Now we shall consider a particle moving around a Taub-NUT attractive center sur- 
rounded by a weak dipolar halo. In Fig. 3 we draw the Poincare section for D = 0.0005, 
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b = 0.14107, a = 1, E = 0.967 and L = 3.75 (direct rotation). We see that the islands of 
stabihty are larger in this case than in perturbed Schwarzschild solution (see Fig. 2); also we 
have new systems of small islands immersed in the chaotic region. We have that the chaotic 
region is smaller in this case than in the perturbed Taub-NUT case. It does look like that 
the presence of the NUT parameter diminishes the effect of the dipolar strength as a chaos 
source. In Fig. 4 we present the section with the same parameters of Fig. 3, except that 
now we have increased the NUT parameter, b = 0.254787. In Fig. 4 we observe that the 
integrable region increases in such a way that the chaotic region has almost disappear in this 
scale. In other words, the NUT parameter has the property of restoring the invariant tori. 
In Fig. 5 we present a magnification of the region around the point (0.0744, 7.37) of Fig. 4, 
we find a small chaotic region around the crossing of the section around the indicated point. 
Although in the present Letter we present results for particular values of the parameters 
involved, we did a rather extended numerical study that supports our conclusions. Figs. 
3, 4, and 5 being representative of this search. The change of sign of the angular particle 
angular momentum does not alter the figures and the change of sign of D introduces only a 
Pr —Pr transformation. 

For the values of the parameters D = 0.0005, b = 0.14107,0.254787 E = 0.967 and 
L = ±3.75, we have that the particles move roughly in the "box" 4.7 < r < 21, —6 < z < 7. 
We take as a measure of error the quantities 

^9a, = m - 9ab)/9:^\, A/ = \{r - f)/n, (§) 

where in these expressions the sum rule of repeated indices does not apply, g^^ and g^i, refer 
to the solutions (0) and (|^), respectively. We find that for the above mentioned range of 
coordinates and the values of parameters used in this Letter the quantities defined in (^ are 
at most of the order of 10~^. Also in this range, the error in the derivatives of the metric 
functions is even smaller (the metric functions are very smooth). We also want to mention 
that the Poincare sections shown in this Letter were computed with an accumulated error 
in the "energy" [cf. Eq. (0)] smaller than 10^^°. 
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In summary the NUT parameter has the rather surprising property of making the mo- 
tion of particles more integrable, it enlarges the region of invariant tori. This is a puzzling 
result since the addition of the NUT parameter and a dipole moment to a black hole makes 
the metric rather complicate, hence one should foresee more complex motion with a greater 
destruction of tori. It happens the opposite. This is again a manifestation of Misner char- 
acterization of the Taub-NUT metric as a 'counterexample to almost anything" . 

The authors thank CNPq and FAPESP for support and to S.R. Ohveira for discussions. 
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FIGURE CAPTIONS 



Fig. 1. Poincare section of test particles moving with angular momentum L — 3.75 and 
energy E = 0.967 in a Taub-NUT geometry with NUT parameter h = 0.28 and (7 = 1. This 
is a typical section of an integrable system. 

Fig. 2. Poincare section for D = 0.0005, E = 0.967, L = 3.75, 6 = 0, and cr = m = 1. The 
two isolated islands around the points (10, 0.05) and (5, 0.075) are parts of the same torus. 

Fig. 3. Poincare section for D = 0.0005, b = 0.14107, a =1 E = 0.967 and L = 3.75. The 
islands of stability are larger in this case than in the geometry without NUT parameter (cf. 
Fig. 2). 

Fig. 4. Poincare section with the same parameters of Fig. 3, except that now b — 0.254787. 
The chaotic region is almost inexistent. 

Fig. 5. The region around the point (0.0744, 7.37) of Fig. 4 is magnified. We find a small 
chaotic region. 
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